Ground state splitting of 8 S' rare earth ions in semiconductors 
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We propose a new mechanism leading to the ground state splitting for the rare earth 8 S ions in 
semiconductor crystals. The resulting splitting is due to three effects, the first is the intra atomic 
4/ — 5d spin-spin interaction, the second one is the spin - orbit interaction for 5d electrons and the 
third one is their hybridization with the valence band states of semiconductors host. The resulting 
splitting significantly depends on the relative position of 5d level with respect to semiconductor host 
band structure. We also discuss different model, already known in the literature, which is also based 
on ion - band states hybridization. For both models, as an example, we present results of numerical 
calculations for rare earth ion in IV- VI semiconductor PbTe. 

PACS numbers: 71.70.Ch, 71.70.Ej 



I. INTRODUCTION 

The present paper is devoted to the theoretical analy- 
sis of the ground state splitting of rare earth (RE) 8 5 ions 
semiconductor crystals. Very well known and studied ex- 
amples of such ions are Eu 2+ and Gd 3+ . It is believed 
that the magnetic properties of these ions are determined 
by the half hlled 4/ shell. According to the Hund's rule, 
for such electron configuration, the ground state should 
be characterized by the total angular momentum L=0 
and the total spin S=7/2. In other words the ground 
state should be 8 S7/ 2 which is eight fold degenerate, in- 
dependently of the crystal environment. Such a model 
explains quite well magnetic susceptibility or magneti- 
zation measurements where the spin of the ion interacts 
with external magnetic field or with the spin of another 
ion and the energies of Zeeman or ion-ion exchange in- 
teractions are of the order of 1 K. 

However, the electron paramagnetic resonance (EPR) ex- 
periments, which probe the system on a much finer en- 
ergy scale, clearly demonstrate that the above picture 
is an approximate one only. It turns out that in crys- 
tals the degeneracy of the ion's ground state is lifted, 
and the nature of the splitting depends on the symmetry 
of the enviroment. For example, for Oh symmetry, the 
case we deal in the present paper with, the ground state 
is split into three levels, doublet, quartet and doublet. 
In accordance with the group theory, the effective spin 
Hamiltonian describing the splitting caused by the ion's 
crystal environment is of the form 
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where the operators equivalents O™ for spin S=7/2 are 
8x8 matrices defined, for example, in Ref. [l. The 
energy levels of Hamiltonian (Q are (— 186 4 — 126 6 )^ 2 \ 
(2b 4 + 166 6 )( 4 \ (146 4 - 20& 6 )( 2 \ where the superscripts 
( 2 ) or ( 4 ) denote the degeneracy of the levels. 
From EPR experiments, if the quality of the samples is 
high, it is possible to obtain not only the absolute values 
of coefficients 64 and be, but also their signs. This may 
be achieved by performing measurements of the same 
sample in different temperatures. The splittings, to- 



FIG. 1: Differences in the splittings of the ground state levels 
for a) Eu 2+ and b) Gd 3+ ions in PbTe crystal. The numbers 
show the degeneracies of the levels. 



gether with levels' degeneracies are shown schematically 
in Fig. n for 64 > and &4 < 0. The energy diagrams in 
Fig. n correspond to the splittings of Eu 2+ and Gd 3+ ions 
in PbTe. This compound is an example we will often re- 
fer in the present paper to, but the calculational methods 
and the general predictions of the theory may be applied 
to any semiconductor containing RE S state ions. From 
the literature we know, that for Eu ion^ &4=129 MHz, 
while for Gd ion 3 64= -110.16 MHz. The coefficient &6 is 
usually about two orders of magnitude smaller and be- 
cause its influence on the ground state splitting is very 
small it will not be analyzed in the following. Notice that 
the signs of 64 for Eu 2+ and Gd 3+ ions are opposite. This 
is rather strange, because at the first sight, except for the 
small relative difference of nuclear charges these two ions 
are very similar and their crystal neighbourhoods are the 
same (six tellurium atoms placed in the vertices of regu- 
lar octahedron). 
The problem of the ground state splitting of rare earth 
8 SV/2 ions is very old. Very comprehensive discussion of 
the possible physical mechanisms leading to the splitting 
is due to Wybourne4i In Ref. using perturbation the- 
ory and results of numerical calculations he analyzed the 
case of Gd 3+ ion in the crystal environment of D^h sym- 
metry. The main idea is, that due to strong spin-orbit 
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interaction for 4/ electrons the higher energy 4/ 7 states 
2S+1 L 7 / 2 with different L and S are mixed into the ion's 
ground state 

| 8 S 7 /2> - s\ 8 S 7/2 ) +p\ & P 7/2 ) + d\ 6 D 7/2 )... (2) 

where s, p, d, ... are numerical coefficients^ Due to 
nonzero p, d, ... the ion is no longer in pure 8 <SV/2 state, 
it may interact with the crystal field and this iteraction 
leads to the ground state splitting. 

The results obtained by Wybourne have been recently 
improved by Smentek et. al£. On the basis of exten- 
sive numerical relativistic calculations for Gd 3+ ion they 
derived the effective spin Hamiltonian for this ion in the 
crystal neighborhood of arbitrary symmetry (relativistic 
crystal field theory). In particular, for Oh symmetry the 
effective spin Hamiltonian may be written in the form of 
Eq. (JTJ where the coefficient 64 reads 

b 4 = |g (^M04) + M13)* (13)4 + &4(15)* (15)4 ) 

The definitions and numerical values of 64 (/ct) 

and X( K *) 4 may be found in Ref. 0. Aq is the crystal 
field coefficient, which, in the simplest six point charge 

model^ is equal Aq = 7 * Jfo where a «0.5 A is the 
atomic length unit, d is the distance between cation and 
anion and e is the electron's charge. Taking Eu and Gd 
in PbTe as example, we see that the relativistic crystal 
field theory alone cannot simultaneously explain ground 
state splitting for both ions because the experimentally 
determined signs of 64 are opposite. It is very improbable 
that replacing Eu ion by Gd ion we change the sign 
of Aq because the neighborhoods of two ions are the 
same. We also do not expect significant changes of the 
parameters t/' 4 ', 64 (/ct) and X^^ 4 because these two 
ions are very similar. 

The inadequacy of a model in which the RE ion interacts 
with the enviroment by electrostatic crystal field poten- 
tial only has been already noticed in the literatur o 5 i 8 i 9 . 
In 1978 Barnes et. al. a , analysing the ground state 
splitting for Gd 3+ ion in different crystals, noticed that 
in insulators there is proportionality between the crystal 
field coefficient Aq and coefficient 64, while in metals 
such proportionality does not exist. They proposed a 
model in which 4/ electrons interact with the band 
states via hybridization processes. Using second order 
perturbation theory with respect to the hybridization 
they constructed effective spin Hamiltonian for the 
ground state of Gd 3+ ion. As the excited states of the 
system they took into account configurations in which 
the number of electrons on 4/ shell changes by ±1, i. e. 
4/ 8 configuration plus one hole in the Fermi sea or 4/ 6 
configuration plus one additional electron in the band. 
According to the Hund's rule, in the excited states 4/ 8 
and 4/ 6 the angular momentum is nonzero. Taking into 
account internal spin - orbit coupling, the authors of 
Ref. obtained effective spin - lattice interaction leading 
to the ground state splitting. 



In Section IV we will rederive this model because we 
think that in its derivation and solution presented in 
Ref. some important points have been missed. 
The main idea of Barnes et. al. is very interesting, 
however the final results strongly depend on the con- 
stants describing the hybridization between 4/ electrons 
and the band states. Due to strong localization of 4/ 
shell, in the literature, it is very often assumed that it is 
very small. For example in calculations of the exchange 
integral between 4/ spin and band carriers in PbEuTc, 
Dietl at. aliiS found that the contribution to the final 
result resulting from the hybridization 4/- band states 
processes is negligible. 

Contrary to 4/ electrons, 5d states of RE ions are very 
extended in space and their hybridization with band 
states is certainly much stronger. In the literature we 
have found examples of successful explanations of mag- 
netic properties of RE ions in semiconductors which are 
based on the assumption that the interaction between 
4/ electrons and band states goes via internal 4/ — 5d 
exchange interaction and 5d - band states hybridization. 
For instance, using such a model, in 1970 Kasuyaii 
explained the Eu-Eu exchange constant in EuO. The 
Kasuya mechanism was used also, as a starting point, 
by Story et. alJti in the theory explaining the Fermi 
energy dependence of the Gd-Gd exchange constant in 
SnGdTe mixed crystals. The 4/ - 5d interaction was 
also invoked by Dietl et. alii^ in calculation of sp — / 
exchange integral between localized Eu spin and the 
band carriers in PbEuTe. 

In the present paper we generalize this mechanism 
by including spin-orbit coupling and crystal field poten- 
tial for 5d states and apply it to the calculation of the 
coefficient 64. In the model, which will be explained 
and discussed in detail in the next sections, the ground 
state splitting occurs due to combined effect of the intra- 
atomic, Heisenberg type, 4/-5d exchange interaction, 
spin-orbit interaction on 5c? orbitals and the hybridiza- 
tion of RE 5d levels with the valence band states. Our 
model takes automatically the spin-orbit interaction in 
the semiconductor host band states. It turns out that the 
functional dependence of 64 on 5c? spin-orbit constant is 
different for semiconductors with strong band spin-orbit 
interaction from those for which this interaction may be 
neglected. For semiconductors with strong band spin- 
orbit effects one may expect that the proposed mecha- 
nism is more effective. One of the important parameters 
of the model is eo, the energy necessary to transfer an 
electron from the valence band to the 5d shell. The mag- 
nitude of the resulting 4/ ground state splitting decreases 
very quickly with eo. That is why the position of 5d level 
with respect to the Fermi energy decides whether the 
mechanism is important or not. For example, as it will 
be discussed in Section III, for Gd in PbTe this energy 
is small, of the order of 0.5 eV— and the resulting 64 is 
of the order of the one observed in experiment, while for 
Eu £q is a few times larger and the calculated splitting is 
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much smaller. 

In the next Section we describe the model and we derive 
approximate, analytical result for 64. This formula en- 
ables us to discuss the salient features of the model, in 
particular the dependencies of &4 on the parameters of 
the theory. In Section III we present details of numeri- 
cal calculations and Section IV is devoted to rederivation 
and discussion of the model by Barnes et. al£. Some 
additional remarks are contained in the last Section. 



II. THE EFFECTIVE SPIN HAMILTONIAN 

Let us consider a semiconductor crystal with one cation 
replaced by RE atom. The unperturbed part of our 
model Hamiltonian describes ground and excited states 
of the system. These two groups of states are connected 
by 5d - band states hybridization, which we treat here 
as a perturbation. Below, we describe these two parts of 
Hamiltonian. 

In the ground state of the system the electrons fill the 
band levels up to the Fermi energy. We assume that the 
band structure is not changed significantly by the pres- 
ence of RE atom. In our model RE atom is treated as 



the host cation atom with additional 4/ and 5c? orbitals. 
There are seven electrons on 4/ orbitals and 5d shell is 
empty. Assuming validity of the Hund's rule, the spin of 
the ion 5 = 7/2 and its angular momentum L = thus 
the ground state of the system is eight fold degenerate 
and is described by M, the projection of spin on a quan- 
tization axis. 

As the excited states we take configurations with one ad- 
ditional electron on the 5d level and one electron less in 
the band. The Hamiltonian describing 4/ 7 5d 1 configura- 
tion of the ion reads: 

H^fbd 1 = -JS ■ s + A t L t ■ s + A e L e ■ s + V cr (4) 

The first term describes the exchange interaction between 
the 4/ spin S and the spin s = \<j of the 5d electron. The 
second and the third terms describe the spin orbit inter- 
action on the bd shell. The operators L t i (i = x, y, z) are 
the angular momentum operators between the 5d states 
of t2 g symmetry and L e i between the tig and ei g states 
respectively. For the definition of tig (dyz, d xz , d xy ) and 
e g (d z 2, d x 2_ y 2) states according to Ref. Qthey have fol- 
lowing form: 
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(6) 



In the case of free ion two spin orbit constants Xt and 
A e are equal, however if the ion is placed into the crystal 
they are in general different £ Finally, the last term in 
the Eq. describes the influence of the crystal field on 
the 5d energy levels and it has the form 



V cr = Dq 



-4 
-4 
0-400 
6 
6 



(7) 



After diagonalizing of the 80x80 matrix H^p^i we ob- 
tain eigenvectors \R) and the corresponding eigenvalues 
cr. The eigenvectors \R) may be expressed in the basis 



\Mcka) 



\R)= ^2 \Mdia){Mdi(x\R) 



(8) 



where —7/2 < M < 7/2 and a — ±i are the projections 
of the 4/ and 5d spins on the quantization axis, respec- 
tively and di are t^g for i = 1,2,3 and e g for i = 4, 5 5d 
orbitals. 

In our model we assume that the 5d levels of RE hybridize 
with the band states. (The hybridization of the 4/ shell 
in the present Section is neglected.) The hybridization 
matrix elements (i?, q\h\M) describe the probability am- 
plitude of a transition from ground state \M) to ex- 
cited state \R, q) with h being one-electron hybridization 
Hamiltonian and the quantum number q = nk describes 
the band state of an electron with the corresponding en- 
ergy e g transferred to bd shell. The wave vector k belongs 
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to the first Brillouin zone and n is additional index neces- 
sary to fully characterize the band state. For bands with 
negligible spin orbit coupling, where the electron spin is 
a good quantum number, n corresponds to the band's 
number and the projection of electron's spin. 
In the second order perturbation theory with respect to 
the hybridization between 5d levels of RE ion and the 
band states, which often is called Schrieffer- Wolff trans- 
formation, we obtain the effective spin Hamiltonian for 
the RE ion in the crystal: 



H 



MM' 



E 

n, q 



(M\h\R, q)(R,q\h\M') 

£0 + (R - tq 



In Eq. © the sum over q runs over all occupied band 
states. The sum in the denominator is the energy of the 
excited state of the system where eo denotes the energy 
necessary to transfer an electron from the Fermi level, 
which we assume to be zero of the energy scale, to the 
lowest energy state of 4/ 7 5d 1 configuration. 
Using Eq. (JHJ) and the fact that 



{Mdiaq\h\M') = 6 M m> (dk<r\h\q) 



(10) 



(9) Eq may rewritten in the form 



h M m> - - zl 

R,q diicri,di2f2 



(d l2 a 2 \h\g) (#|dqai) {Mdga^R) (R\M'd l2 a 2 ) 

(Q + tR - tq 



(11) 



r 



Eq. Ijlljl is the main result of the present paper. If we 
know the values of matrix elements Hmm 1 then, compar- 
ing them to Eq. QJ, we obtain the coefficient 64. In the 
next section we present details and results of numerical 
calculations for the example case of PbTe semiconductor, 
here we derive approximate, but analytical formula for 64 
valid to the fourth order with respect to the intra-atomic 
spin-orbit coupling. This analytical formula will enable 
us to discuss and understand the general properties and 
dependencies of 64 on different parameters of the model. 
Let us take an arbitrary state characterized by wave 
vector k from the first Brillouin zone and let us de- 
note by {k } the set of states which may be obtained 
from ko by symmetry transformations of Oh, including 
Kramers conjugation. More precisely, we take an arbi- 
trary wave vector ko from the first Brillouin zone, we 
find all the wave vectors which may be obtained from 
ko by symmetry operations and the set {ko} contains 
the Kramers conjugate pairs corresponding to these vec- 
tors. On the symmetry ground, the matrix Zd il<Jl ,d i2 <T2 = 
S 9 e{k }( c ^i CT il%)( l 3 , H^2C''2), appearing in Eq. (JTTJ, 
must have following form: 

ZdiKri,da<ra = Vth + Peh + Qt^t ' CT + 9eL e • CT (12) 

where the matrices /( and I e are unit operators in t 2g 
and e g subspaces, respectively, i. e. : 



h = 



1 

10 

10 







h = 









1 

1 



(13) 



The coefficients pt,p e ,qt,qe depend on the band structure 
and on the quantum number ko. If we neglect the band 
spin-orbit interaction than q t and q e disappear. However, 
in many cases, for example for PbTe, such an assumption 



is unjustified and, as we will see below, q t and q e , in some 
sense, play the role of ion's spin-orbit constants. 
Let us notice that e q = ek =const for all q G {ko}. 
Than the effective Hamiltonian Hmm 1 , Eq. pif) . may be 
rewritten as the trace over di and <j degrees of freedom 
of the product of two matrices: 



Hmw =^J2Tr di(T ZQ{M,M') 



(14) 



where the matrix Q(M, M') is defined as 
n , M „/, \- (MdiKnlRHRlM'djM) 

Qd ll a 1 ,d i2 a 2 \ M , M ) — — > ■ 

^ eo + e R - e ko 

(15) 

and the factor 1/48 is necessary to take into account mul- 
tiple counting of states (the Oh group has 48 elements). 
Eq. H15|) may be expanded in the series of At and A e . De- 
noting H\ = A 4 Lt • s + A e L e • s the expansion is according 
to the formula: 



1 



E - 



7 JS ■ (T + V cr + H\ 



= G a J2 (~HxGoT (16) 



n=0 



where E = eo — e q + 7/4 J + 4Dq. Notice that the low- 
est eigenvalue of the operator — ^JS • er + V cr equals 
— (7/4 J + ADq) and the subtraction of this term is in 
accordance with the definition of eo. The operator Go = 
1/ (E+ |JS • er + V cr ) is equal: 



Go = {A t + B t S • cr) I t + (A e + B e S ■ a) I e 



(17) 



with 



eo 



J 



(eo - e q ) (e - e q + 4J) 



(18) 
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2 (eo - e q ) (e - e g + 4 J) 

6 - e 9 + 10J»g + f J 

(e - e q + WDq) (e - e g + 10D« + 4J) 

(19) 



e 2 (e - e q + WDq) (e - e q + IQDq + 4 J) 

Calculations of traces in Eq 1)14)1 is a simple but very te- 
dious task, which may be simplified by use of a computer 
program. We collect terms proportional to 

S x + St +Si = ± (0° + 50j) + ^ (20) 

and we obtain the coefficient 64, up to the fourth power 
in ion's spin-orbit coupling: 

bi = ~ B t B ^ {qtBtXtXe + q e (B t X 2 t - 2B e \ 2 e ) 

ko 

-^A e [2qt (A t - A e ) B t X\ + 8q e (A e B t - A t B e ) X t X e 

+ Pt A t (3B t X 2 t - 2B e X 2 e ) + p e A e (B t X 2 - 2B e X 2 e )] } ■ 

(21) 

Let us discuss the main features of the obtained formula. 
First, let us notice that b^~Q for A e =0. This is the gen- 
eral result valid for all orders of perturbation theory. It is 
connected with the fact that if A e =0 there is no spin orbit 
coupling between i 2g and e g orbitals and this coupling is 
the only one which mixes these two groups of states in 
Hamiltonian H\. The operator Go, Eq (|17fl does not 
mix tig and e g states too and, consequently, the same 
holds for operator Qd il<Tli d i2<T2 (M, M'), see Eqs. 1)151161 . 
Notice, that for A e =0, in Eq. ltPo|) . in the subspace of 
e g states, only the zeroth order term survives and this 
term will certainly not lead to terms in the effective spin 
Hamiltonian which have Oh symmetry. Moreover, due 
to this decoupling, the trace in Eq. 1)14)1 does not depend 
on q e . Thus, concerning the orbital degrees of freedom, 
we may limit the considerations to the ti g subspace. But 
in this subspace the problem is completely symmetrical 
with respect to operations of full rotational group be- 
cause L tx -,Lty,Lt z satisfy the angular momentum commu- 
tation relations for L=l and there is no reason to obtain 
Oh terms in the effective Hamiltonian. We conclude that 
the nonzero spin-orbit coupling A e connecting the ti g and 
e g states is the most important parameter of the model. 
Secondly, notice that the order of ion's spin-orbit cou- 
pling, A, in which we obtain nonzero 64 depends on the 
band spin-orbit coupling. If the spin-orbit coupling is ab- 
sent in the band, q t = q e =Q, the lowest order terms are 
proportional to A 4 . For nonzero band spin-orbit coupling 
they are proportional to A 3 . In this sense the band spin- 
orbit coupling plays the role of ion's spin-orbit coupling. 
The formulas for A t , B t , A e , B e show that the coefficient 



&4 very quickly decays with the excitation energy eo- On 
the other hand, these formulas suggest also that the most 
important contribution to the final results comes from the 
hybridization of 5d level with the states from the vicinity 
of the Fermi level. 

Finally let us notice that the crystal field potential for 
Dq > diminishes 64. This may be understood on the 
basis of the preceding discussion devoted to the role of A e . 
Enlarging the energy distance between t2 g and e g states 
leads to the effective decreasing of the coupling between 
these two groups of states. 



III. RESULTS 

In this Section we show the example numerical calcu- 
lations for RE ion in PbTe semiconductor. 
The electron band wave functions and band energies e q 
are calculated according to the tight binding model de- 
veloped in Ref. Q. In this model the band states are 
build from p and s orbitals of Pb and Te. For a given 
momentum k belonging to the first Brillouin zone, the 
tight binding Hamiltonian is diagonalized in the basis of 
16 functions of the form 

tf£(r) = -j= £ e**"*tf'(T - R c/a ) \a), (22) 

where tp c J a (r — R c / a ) with i — p x ,p y ,p z , s are the cation 
or anion atomic orbitals centered on the lattice sites R c 
or R a , respectively, iV c is the number of cation sites, and 
\a) with a — ±i is the two dimensional spinor. After 
diagonalizing the Hamiltonian matrix, for a given k we 
obtain the band energies e„k and corresponding eigen- 
functions 

\q) ee |nk> £ (23) 

i a p—c,a 

where index n = 1, ...16 labels the band number. 

The knowledge of the amplitudes a£ icr enables us to cal- 
culate the hybridization matrix elements (dia\h\q). We 
assume that there is only an overlap of 5d RE orbitals 
with six neighboring anions. The necessary values of in- 
ter atomic matrix elements (dAhltfj) are calculated ac- 
cording to Table 20-1 of Ref. [jja and can be expressed 
usin g th ree constants V s da, V p da, and Vpdrr defined as in 
Ref.]!! 

h 2 r 3 / 2 

Vldm — Vldm (24) 

mod'' 

where rjsdcr = -3.13, rj pda = -2.95 and r\ v ^ = 1.36, m 
is the bare electron's mass, d = 3.2 A is the cation-anion 
distance in PbTe and is a fitting parameter related to 
the radius of 5c? RE orbital and it is of the order of 1 A. In 
calculations we put rd=2.5 A and we obtain V s dcr=-1.62 
eV, V^,d CT =-1.51 eV and V^7i-=0.70 eV. These values are 
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200 




Spin-orbit coupling X ( (eV) 

FIG. 2: Dependence of 64 on the spin - orbit constant At 
for A e =0.1 eV, J=0.25 eV, e =0.35 eV and for two different 
values of crystal field parameter WDq. 

close to the ones used in Ref. ^| in calculations of EuTe 
band structured 

The summation over ko in Eq. 1)11(1 is replaced by the 
integration over the Brillouin zone according to the for- 
mula: 

ko Jbz (2tt) 

where V is the volume of the crystal. 
The value of eo, the energy necessary to transfer an elec- 
tron from the top of the valence band to the 5<i level of 
RE ion is one of the most important parameters of the 
theory. According to Ref. the 5d level of Gd in PbTc 
lies about 0.2 eV above the bottom of conduction band. 
Adding the value of the energy gap in PbTe which is equal 
to 0.2 eV we obtain e =0.4 eV for Gd in PbTe. From the 
resonant photo-emission spectroscopy experiments^ we 
know that 4/ level of Eu in PbTe lies near the top of 
the valence band while for Gd it is placed about 10 eV 
below. That is why we expect that eo for Eu is larger 
than for Gd. The precise value is not important because 
as we will see the value of the coefficient 64 decays very 
quickly with increasing eo and if it is bigger than 1 eV the 
contribution to 64 from the present mechanism becomes 
negligible. 

In the literature we haven't found the values for the spin 
orbit constants At and A e for RE ions in PbTe. One may 
find only value of single spin orbit constant A. The exist- 
ing data are not very consistent, however. In Ref. 19 we 
find A=0.08 eV for Eu and A=0.13 eV for Gd. According 
to more recent theoretical relativistic calculations^, 5d 
spin orbit constant for free lanthanide ions is of the order 
of 0.06 eV. In our calculations we take A e =0.1 eV and in 
Fig. [21 we present the results for b\ as a function of X t 
for eo=0.35 eV, J=0.25 eV and for two different values 
of lQDq = and 10L>q=0.5 eV. We clearly see that the 
nonzero crystal field splitting of 5d strongly decreases 64. 
Taking into account the theoretical calculations 21 of A t 




FIG. 3: Dependence of 64 on the transfer energy eo for 
A e =0.1 eV, A e =0.05 eV, J=0.25 eV and for two different val- 
ues of crystal field parameter WDq. 

and A e performed for 3d ions in crystals of NaCl crystal 
structure, we expect that also in the present case the real 
At should be less than A e . 

In Fig. |2| we show the decay of 64 with increasing of eo for 
At=0.05 eV and for two values of crystal field parameter 
Dq. 

As it is clear from the above discussion the main prob- 
lem in calculations is the lack of knowledge of precise 
values of number of necessary parameters of the model. 
However, it seems that the estimations made above sug- 
gest that the proposed mechanism gives the ground state 
splitting of the right order of magnitude and, moreover, 
in the case of Gd in PbTe it gives proper sign of the co- 
efficient 64. Also, from Fig. we see that it is enough 
to increase eo less than 1 eV to significantly decrease 64, 
which although negative, becomes very small and may 
be neglected. It means that for such cases, for example 
for Eu in PbTe, the main contributions to the splitting 
are due to other mechanisms, like relativistic crystal field 
theory which generates positive sign of 64. 
The above considerations indicate that the proposed 
model should be taken into account in theoretical analy- 
sis of EPR spectra. 



IV. BARNES, BABERSCHKE AND 
HARDIMAN MODEL 

As it has been discussed in the Introduction the main 
idea of the model proposed by Barnes et. &\~ is to con- 
sider the excited states of the system in which the number 
of electrons on the ion's 4/ shell changes by ±1. Let us 
concentrate in this Section on processes which lead to 
4/ 7 <-> 4/ 6 transitions. The ground state of the system, 
I S S(J Z )), is the same as in Section II, namely the RE 
ion in 4/ 7 configuration plus the Fermi sea of electrons. 
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This eight fold degenerate state of the system is char- 
acterized by -7/2 < M < 7/2 - the projection of 4/ 7 
spin 7/2 on a quantization axis which we take along the 
(001) crystallographic direction. In the excited states we 
have the ion in 4/ 6 configuration plus one additional elec- 
tron above the Fermi level characterized, like in Section 
II, by the quantum number q. Assuming the validity of 
the Hund's rule for 4/ 6 configuration, L=3 and S=3, the 
Hamiltonian for the ion in excited state reads: 

H = (A 4 L 4 + A5L5 + A25L25 + A45L45) ■ S + V cr (26) 

In the above equation the indices 2,4,5 correspond to the 
decomposition of the D 3 representation into irreducible 
representation T 2 , I\4 and T 5 of the cubic groupi For 
example L4 is the angular momentum operator between 
base functions of T4 representation and L45 between T4 
and T$ base functions. S is the spin operator of the length 
S=3 and V cr is diagonal matrix describing the crystal 
field potential. 

Let us stress that from the physical point of view this is 
rather unreasonable to describe spin - orbit interaction 
using four different constants A4, A 5 , A25 and A45. The 
4/ shell is very localized and contrary to 5d orbitals the 
influence of the neighbourhood should be negligible. The 
only reason of introducing four constants instead of a sin- 
gle one is, that it is easier to understand the influence of 
the crystal field on the coefficient 64. 
The important difference between our formulation of the 
model and that of Barnes et. al. 9 is in the form of hy- 
bridization elements. We propose that 

(L z S z q\H\ M) = (27) 
E ]J 7/2+ 7 2aM 5s»,M-Mh\<t>-^). 

cr=±\ 

The state \L z S z q) is the excited state of the system in 
which the projection on the quantization axis of the to- 
tal angular momentum and spin of the ion are L z and 
S z , respectively, and there is one additional electron 
characterized by q above the Fermi energy. The ele- 
ment {q\h\<p—L z a-) describes hybridization between band 
state q and the 4/ spin orbital 4>^L z a- The coefficient 

(-l) L * +1 y^ 7/2+ 7 2CTM , omitted in Ref. 9, may be derived 
using explicit forms of antisymmetric many electron func- 
tions for ion's states \L Z S Z ) and | M)2£. 
The further steps of calculation of the effective spin 
Hamiltonian are, with minor modifications, very similar 
to those described in the previous sections. Because of 
much higher complexity of the angular momentum alge- 
bra for L=S=3 it is not easy to derive formulas analogous 
to those from Section II. That is why the conclusions are 
based on numerical calculations only, performed for PbTe 
crystal, the example semiconductor we use in the present 
paper. 

In general, the results are similar to those for the model 
described in Section II. First, we observe decrease of mod- 



ulus of 64 with the increase of eo, see Fig.0] In this Sec- 
tion eo denotes the energy necessary to transfer an elec- 
tron from 4/ shell to the Fermi energy level. Secondly, 
for eo an d other parameters of the model kept constant, 
64 changes with the crystal field potential parameter A cr , 
however as we see in the insert of Fig.0]these changes are 
not very fast. A cr is defined here as the total splitting of 
4/ 6 manifold in the crystal field, i. e. A cr = E-p 4 — £r 2 • 
The mechanism of changes of 64 with A cr is similar to the 
one analysed in Section II. From the performed numerical 
analysis it turns out that 64 =0 for A45=0. The relative 
energy positions of T2, T4 and T§ states, which change 
with A cr are important, because they decide about the 
effectiveness of the transitions between them caused by 
spin orbit interaction. However, comparing to 5d, in the 
case of 4/ shell this influence is rather small. This is 
connected with completely different ratio of spin-orbit to 
the crystal field strength. For the 4/ orbitals the effects 
caused by the crystal field are much smaller than those 
due to spin orbit interaction. 

This main difference between our results and those by 
Barnes et. al. is that the authors of Ref. overlooked 
the fact that even for zero crystal field potential we ob- 
tain nonzero ground state splitting of the RE ion. This 
splitting, however, may be achieved only if the symmetry 
of the band wave functions and, consequently, the proper 
symmetry of hybridization matrix elements is taken into 
account. If these symmetries are neglected and the 4/ 
- band hybridization is described by a single constant 
Vfk, (see Ref. then, apart from the crystal field po- 
tential, the rest of the system's Hamiltonian is invari- 
ant with respect to operations of rotation group and the 
ground state remains degenerate. Then, indeed, in the 
first approximation, the splitting will be proportional to 
the crystal field and this is the only contribution calcu- 
lated in Ref. 0. 

On the basis of our analysis, see Fig. 0] it is clear that 
in the case of 4/ - band states hybridization, the most 
important parameter of the model is not A cr but eo, the 
energy necessary to transfer an electron from 4/ shell to 
the conduction band. 



V. CONCLUSIONS 

In this paper we have proposed the model leading to 
the splitting of the ground state of S S rare earth ions in 
crystals. The main ingredients of the model are 4/ — 5d 
exchange intereaction, spin orbit coupling for 5d electrons 
and their hybridization with the band states. The numer- 
ical calculations have been performed for PbTe, a semi- 
conductor from IV- VI group of compounds. We have also 
limited considerations to the Oh symmetry of the ion's 
neighbourhood. Of course, the model may be applied to 
any semiconductor and, with slight modifications in hy- 
bridization matrix elements, to lower symmetry cases. 
In Section IV we have discussed the model proposed by 
Barnes et. al. In particular, we have shown the impor- 
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bridization matrix. If the symmetry of the problem of 
this kind ist lost due to too big simplifications in the 
model's formulation we may also loose important, from 
the physical point of view, class of solutions, as it has 
happened in Ref . |9|. 

The general conclusion resulting from the analysis of both 
models is inaccordance with the point of view presented 
in Ref. namely that the ground state splitting of 8 S 
ions in crystals is not governed by their internal proper- 
ties only, but the position of the ion's 4/ and 5d levels 
relative to the host crystal band structure is at least of 
equal importance. 

The model introduced in Section II, applied to Gd 3+ and 
Eu 2+ ions in PbTe explains origin of opposite signs of the 
coefficients 64 observed in experiment. 



FIG. 4: Dependence of &4 on the transfer energy eo for 
A 4/ =0.0276 eV, V sdtT =-0.2 eV, V pdlT =-0.2 eV and 1^=0.1 
eV and the crystal field splitting A cr =0. The insert shows 
the dependence of 64 on A cr for eo=0.5 eV. 



tance of the careful treatment of symmetry of the hy- 
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